Orbital ac spin-Hall effect in the hopping regime 
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The Rashba and Dresselhaus spin-orbit interactions are both shown to yield the low temperature 
spin-Hall effect for strongly localized electrons coupled to phonons. A frequency-dependent electric 
field E(w) generates a spin-polarization current, normal to E, due to interference of hopping paths. 
At zero temperature the corresponding spin-Hall conductivity is real and is proportional to u> 2 . At 
non-zero temperatures the coupling to the phonons yields an imaginary term proportional to u. The 
interference also yields persistent spin currents at thermal equilibrium, at E = 0. The contributions 
from the Dresselhaus and Rashba interactions to the interference oppose each other. 

PACS numbers: 72.25.Dc, 73.21.-b, 71.70.Ej 



Introduction. The possibility to control electron spins 
by an electric field, due to spin-orbit (SO) interactions, 
has obvious potential applications. Present-day research 
has focused on the spin-Hall effect: An in-plane electric 
field applied on a two-dimensional electron gas creates 
an in-plane spin-current, flowing in the perpendicular di- 
rection. In a two-dimensional electron gas confined to 
the x — y plane, the spin-current is the flux of electrons 
with spins polarized along z. The spin-Hall conductiv- 
ity, defined as the ratio between the spin-current and 
the electric field, takes a universal value [1] in a pure, 
homogeneous infinite electron gas, but disappears [2] in 
the presence of arbitrarily small static disorder. Conse- 
quently, the spin-Hall effect recently observed in a two- 
dimensional hole system [3] is apparently due to the sam- 
ple's edges. 

Spin-orbit interactions in two-dimensional electronic 
systems originate from bulk inversion asymmetry (due to 
the Dresselhaus [4] term), or from the structural inversion 
asymmetry of the potential confining the electrons to the 
plane (the Rashba [5] term). While many of the theo- 
retical studies concentrate on either term [6] , it has been 
noted [7-9] that there is, in fact, a competition between 
the two terms, implying the intriguing experimental pos- 
sibility to control the direction of the spin-polarization 
flow by modifying, e.g., the asymmetry of the confining 
potential. 

Whereas the spin-Hall effect of itinerant electrons in 
the diffusive regime seems to be well understood, far 
less has been done concerning its realization in strongly- 
localized electronic systems coupled to a phonon bath 
[10]. To lowest order, the SO interaction appears as a 
2x2 phase- factor matrix multiplying the hopping ampli- 
tudes [11,12], 



Here, Ju — Ju is the overlap of two wave functions local- 
ized at sites i and i (which can be chosen to be real) , er is 
the vector of Pauli matrices and the vector du = —da is 
calculated below for the combined Rashba and Dressel- 
haus interactions. Therefore, the spin-Hall effect in insu- 
lators is due to interference of hopping paths, reminiscing 
the origin of the ordinary Hall effect in insulators [13]. 
This interference also leads, as we show below, to the ap- 
pearance of 'persistent spin-currents', flowing at thermal 
equilibrium in the absence of any external fields, and sur- 
viving the coupling to the phonon bath - again in anal- 
ogy with the ordinary Aharonov-Bohm persistent charge 
current in the hopping regime [14]. However, this spin- 
persistent current does not cause any local spin accumu- 
lation [15]. To achieve the latter accumulation, and to 
establish a spin-Hall polarization, a frequency-dependent 
electric field, E(w), is required. Below we calculate the 
spin-Hall conductivity matrix of a three-site triad (the 
smallest cluster of localized electronic sites in which the 
SO interaction is manifested), and show that it consists 
of off-diagonal matrix elements alone, such that the spin- 
current density in the plane, j p (w), always flows normally 
to the field, 
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At zero temperature and low u>, we find that the spin- 
Hall conductivity, cr^ all , is real and quadratic in u>, and 
is robust against the self-averaging effect of our highly- 
disordered system. As the temperature is increased, the 
coupling to the thermal bath comes into play, yielding 
an imaginary term in <7g all which results from dissipa- 
tion and is linear in lo. 

As is the case for itinerant electrons [7-9] , the Rashba 
and the Dresselhaus terms work one against the other: 



1 



The interference phase induced by the SO interaction is 
found to be proportional to the difference between the 
two respective contributions. This opens the possibility 
to reverse the sense of the equilibrium persistent spin- 
current, as well as the direction of the spin polarization 
induced by an external electric field, in the strongly lo- 
calized regime. 

Spin-orbit interaction in the hopping regime. The com- 
bined Dresselhaus and Rashba SO interactions in a two- 
dimensional electron gas may be written in the form 

Hso = U p • a, (3) 

where U p = [a D p x + a R p y , -{a R p x + a D p y )] is a two- 
dimensional vector, consisting of the in-plane momentum 
components, and the SO coefficients an (for the Dressel- 
haus term) and cur (for the Rashba one). To lowest order 
in the coefficients ctD,R, the hopping matrix elements be- 
tween two states localized around sites i and I are given 
by Eq. (1), with 

d ti = [ PD R X U + p R R y u , ~{prR x u + PdR v u )1 (4) 

where pd,r = mctc^ (m is the electronic mass, and h is 
taken as 1) and = Ri — R,£ is the difference of the 
radius- vectors of sites i and £ [16]. 

Transport properties in the hopping regime are cus- 
tomarily obtained as an expansion in the hopping ampli- 
tudes, since the strong localization regime is character- 
ized [17] by \J U \ < \e u \, where e u =e i -e t , e i being the 
single-particle energies of the localized states (assumed to 
be randomly distributed). The single-bond conductance 
(oc Tr{ JuJu}, where the trace is in spin space) is then 
independent of the SO coupling. Spin-orbit interactions 
affect the hopping at third order in the J's (and beyond). 
At third order, one needs to consider the interference be- 
tween the direct hopping path i — ► t and the indirect one, 
i — > m — > £. To second order in the SO coefficients, the 
latter path involves 

- i[da + di m x A mi ] ■ tr). (5) 

At third order in the J's, the interference contribution to 
the conductance comes from traces like Tr{ Jn Ji m J m ,}. 
However, simple algebra shows that, to second order in 
this trace is independent of the SO terms. In con- 
trast, the interference terms do produce a current of spin 
polarization, which results from the SO interactions and 
is quadratic in the d's. As discussed below, this spin- 
polarization current requires the calculation of traces like 

J Urn = -iT?{crJiiJemJmi}- (6) 
The cornerstone of all our results is the observation that 

J Urn = 4:JuJtmJmi(PD ~ (?) 



where the vector A. = R„« x R^/2 is along the z direc- 
tion, and | .4.| is equal to the area spanned by the three 
sites, i, I and m, in the x — y plane. It is also seen that 
the two SO interactions considered here have opposite 
effects on the interference phase. 

Spin polarization and spin-polarization current density 
in the hopping regime. In our calculations we employ the 
usual electron-phonon Hamiltonian pertaining to elec- 
trons in the hopping regime, H = Ho + V, in which 

U Q=Y1 E n n na + E <^& q & q , = c Lc„ CT , (8) 

tig q 

where c) n(J (c nr7 ) creates (destroys) an electron of spin in- 
dex a on the state localized at site n, u q is the phonon 
frequency, and 6 q , 6 q are the phonon operators. The 
Hamiltonian V includes the hopping, the coupling to the 
phonons [14], and the external electric field, 

^ = ^""^ ^""^ 4n' Qnn' e ^ ^ c na c n' a' ! (9) 

nn' (7(7* 

where the spin-dependent hopping is defined in Eq. (1), 




yields the effect of an on-site electron-phonon cou- 
pling (Uq™ = — t)q , with representing the 
electron-phonon coupling on site n) on the hopping, and 
$ nn , is due to an in-plane external electric field, with 
iw$ nn ,(o;) = eE-R nn / [18]. 

The spin polarization operator at site n is given by 
P n = Tr{a z n n<J }. The spin-polarization current density, 
j p , is then defined in analogy with the current density 
due to an electrical dipole-moment [13], 

j p = (l/S)J2(dPn/dt)K n , (11) 

n 

where S denotes the area of the system. Below we con- 
centrate on j p of a single triad, and therefore this area is 
replaced by \ A\. The average, (...}, is calculated with the 
Hamiltonian H assuming that each site is in contact with 
a grand-canonical electron reservoir [13,14]. The tempo- 
ral derivative dP n /dt is found by calculating i\H,n na ] 
up to first-order in the electric field. This yields for the 
average (carried out with the Hamiltonian H) 

(dP n /dt) = - E < 7 ™'' - 
n' era-' 

-i*nn'(*)«'+tfn)>. (12) 

with 

4n' ~ i-Q nn' c na J nn' c n' cr' ■ (13) 
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We calculate (dP n /dt) to third-order in the hopping, and 
up to first-order in the electric field. Detailed calcula- 
tions [19] show that + to third-order in the 
hopping, vanishes. The entire contributions to the polar- 
ization rate comes from — 1%,*) in (12), which will 
be found up to first-order in the electric field. 

Persistent spin currents. The thermal average — 
In'n) i s non-zero even when the electric field is absent. 
This is a consequence of the trace Eq. (6), which includes 
<x: Since the SO interaction conserves time-reversal sym- 
metry, it cannot produce persistent charge currents, but 
it does lead to persistent spin-current. Indeed, ignor- 
ing the electron-phonon coupling, and using third-order 
perturbation theory in the hopping matrix elements, the 
contribution of the bond i — I in our triad (at E = 0) 
becomes 

E *^M° - = U&m E 7^7-. ( 14 ) 



per 



where ft is the Fermi function of the site i occupation, 
and X) P cr stands for the sum over the three permutations 
i,£,m — > £,m,i — > to, i, £. 

Hence, there is a persistent spin-current flowing around 
the triad at equilibrium, even at zero temperature, with 
its sense being determined by the relative locations of 
the site energies compared to the Fermi level. This find- 
ing is in a complete analogy with the persistent charge 
current flowing in response to an Aharonov-Bohm phase 
[14]. The analogy continues when the coupling to the 
phonon bath is switched-on: That coupling induces a 
Debye- Waller factor multiplying the result (14), and also 
a 'counter' spin-persistent current, which flows in the op- 
posite direction to the current (14), and vanishes at zero 
temperature (for details, see Ref. [14]). This equilibrium 
persistent current does not lead to spin accumulation: 
The current flowing from i to £ is minus the one flowing 
from i to to. [When £ is interchanged with to, the sign 
of A. is reversed, see Eq. (7).] Consequently (dPi/dt)\ cq 
vanishes (an analogous calculation shows that at equilib- 
rium, (Pi) |oq vanishes as well). Equilibrium spin-currents 
which do not lead to spin accumulation exist also in con- 
ductors lacking inversion symmetry [15]. The question 
whether they are amenable to an experimental detection 
(like the charge persistent currents) in small coherent 
mesoscopic structures is left open. 

Spin-Hall effect in the hopping regime. When an 
in-plane electric field is applied, the polarization rate 
(dPi/dt) is non-zero. We first consider it for our three- 
site cluster in the absence of the coupling to the phonon 
bath, employing linear response theory (with respect to 
the field) and third-order perturbation theory (with re- 
spect to the hopping amplitudes). The result is 



(dPi/dt){uj) 



J (7(7 



2J< 



i£m 



E 

per 



fi - ft 



(15) 



Namely, the average polarization rate at site i is driven 
by the potential difference across the bond £ — to. This 
implies that inserting Eq. (15) into Eq. (11) produces 
Eq. (2), with the spin-Hall conductivity, cr^ all , given by 



okuiM = WeLo 2 (p 2 R -p 2 D )\A\Tt im (uj) 



where 



■ Urn 



JilJimJmi \ A 



~i£ e £m e mi per C U 



(16) 



(17) 



In deriving this result, we have discarded resonance tran- 
sitions in which the frequency u) compensates the site en- 
ergy differences. Otherwise, the sum in Eq. (17) would 
have been augmented by terms of the type fe)5(e ie ± 
We will find below that the electron-phonon cou- 
pling gives rise to an imaginary part in the spin-Hall con- 
ductivity, that is of the same order in ui, and which origi- 
nates from energy-conserving delta-functions of the type 
S(e ie ± uj q ± lo). Since the latter are more likely to be 
satisfied for the randomly-distributed site energies, the 
present discussion is confined to the result (16). 

The result (2) demonstrates in a nutshell the origin of 
the present-day interest in the spin-Hall effect. While an 
electric field along, say the x direction, will drive a charge 
current in the same direction, the spin-current density in 
this case will be along the y direction, leading to a 'sep- 
aration' of spin and charge. 

The spin-Hall conductivity (16) vanishes at small fre- 
quencies and tends to zero as lo~ 2 at very high ones [since 
when e 2 ^ <C to 2 the sum in Eq. (17) vanishes]. Thus, a 
dc electric field is incapable of producing the spin-Hall 
effect. This is in accordance with the behavior found in 
itinerant electron systems [2]. This conductivity, which 
depends upon the temperature through the Fermi func- 
tions, remains finite at zero temperature when two of the 
three site energies are below or above the Fermi energy. 
Then, at frequencies smaller than the typical site energy 
differences, and assuming that all the hopping matrix el- 
ements Jn have the same sign [12], T c Um takes a definite 
sign, independent of the locations of the site energies rel- 
ative to the Fermi level. Consequently, the self-averaging 
effect of the macroscopic hopping system, which involves 
many triads will not wash out cr^ all . 

The effect of the coupling to the phonons on the 
above results may be divided into three. Firstly, there 
is the overall Debye- Waller exponent, due to loss of co- 
herence. This factor exists even at zero temperature. 
(For brevity, it will not be presented below.) Secondly, 
virtual electron-phonon processes contribute additional 
terms to r* (m , which are proportional to the electron- 
phonon coupling and involve complicated combinations 
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of site Fermi functions and Bose occupation numbers of 
the phonons. We do not present these terms since in 
the weak electron-phonon coupling they are smaller than 
the purely electronic ones given in Eq. (17). Thirdly 
and most importantly, real (energy-conserving) clcctron- 
phonon transitions give rise to an imaginary part in the 
Hall-conductivity, which now becomes 

^nH = 16e(p 2 fl - p 2 D )\A\(^T<r £m (co) *r£>)) (18) 

with 



sinh(/3w) 



nor 



G%1 Ji m. J rr 

e 2 Jm 



(19) 



where Ga oc exp[— /3(|e i£ | + \e { \ + \e e \)/2] is the usual 
temperature-dependent hopping conductance of the i — I 
bond, and (3 is the inverse temperature. The result (19) 
is derived in the small uo limit, assuming the frequency to 
be smaller than the site energy differences. Hence, r op 
vanishes at zero temperature. 

To obtain the temperature dependence of the Hall- 
conductivity (18) at low frequencies, we consider the sit- 
uation in which the leading electrical conductance of our 
triad takes place along the bond i — £, and site m provides 
the interference path necessary for the spin-Hall effect. 
Thus we imagine e i and to be below and above the 
Fermi level, but close to it, while e m lies far away from 



the Fermi energy. In that case 



J it J tvn, Jmi 1 1 ^ \ 



while IX 

Jit = 



t\^mi 



r «m7«sinh(/3w)|e^|exp(-/3|e l£ |), where 
2 7V(|e i£ |) is the density of phonon states at 
|e^| multiplied by the electron-phonon matrix clement 
squared. It then turns out that one may define a charac- 
teristic frequency of the system [18], 



MT)=JuP\e u \e-^, 

which vanishes at zero temperature. 
spin-Hall conductivity is then 

^HallM XU[U-WU{T)]. 



(20) 

For Puj < 1 the 



(21) 



At very low temperature, this conductivity is propor- 
tional to lo 2 ; As the temperature is increased, real 
electron-phonon processes give rise to an imaginary part 
due to dissipation, which is linear in uj. 

In conclusion, we have shown that the competition be- 
tween the Rashba and the Dresselhaus SO interactions 
opens the possibility to control the magnitude and the 
phase (compared to the driving ac field) of the spin-Hall 
current for localized electrons in the hopping regime. 
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